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Abstract 


Yield  strength  obtained  fiom  quasi-static  strength  data  for  rolled  homogeneous  armor 
(RHA)  was  combined  with  dynamic  strength  data  for  2-in  (51 -mm)  RHA  to  generate 
Johnson-Cook  parametCTs  for  2-in  RHA.  One  parameter  was  fixed  based  on  the  quasi¬ 
static  strength  data,  and  a  least-squares  method  was  used  to  fit  the  others  individually. 
The  fit  was  tested  with  CTH  by  simulating  the  penetration  of  stacks  of  2.5-in-thick 
(63.5-mm)  RHA  plates  (the  closest  available  experimental  data).  Parameter  analysis  and 
comparison  of  the  simulations  to  experiment  substantiated  the  approach. 
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1.  Introduction 


Class  1  rolled  homogeneous  armor  (RHA),  designed  for  maximum  penetration  resistance,  is 
available  in  thicknesses  from  1/4  in  (6.35  mm)  up  to  6  in  (152.4  mm).  Military  specification 
MIL-A-12560H  (U.S.  Department  of  Defense  1991)  allows  a  wide  variation  in  hardness  over  the 
range  of  available  thicknesses  as  well  as  within  each  thickness  group  (Figure  1).  Since  hardness 
is  an  indicator  of  several  material  strength  properties,  a  significant  variation  in  material 
properties  exists  over  the  range  of  thicknesses  of  available  RHA  and  to  a  lesser  extent,  within 
each  thickness  group. 
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Figure  1.  RHA  Hardness  Variations  Specified  by  1VHL-A-12560H. 


To  properly  model  the  ballistic  performance  of  RHA,  consideration  must  be  given  to  these 
propeily  variations.  The  wide  variation  in  material  properties  across  the  spectrum  of  available 
thicknesses  suggests  that  each  thickness  should  be  separately  evaluated  to  obtain  valid  strength 
parameters.  Further  complications  exist  (e.g.,  manufacturing  lots  and  through-the-thickness 
hardness  variations).  However,  these  eomplexities  are  avoided  in  the  present  work  by  assuming 
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that  the  variations  in  properties  for  a  particular  thickness,  as  allowed  by  the  thickness  group,  are 
negligible.  That  is,  for  an  RHA  plate  that  confoims  to  MIL-A-12560H,  specifying  its  thickness 
is  sufficient  in  identifying  its  properties. 


The  shock  physics  code  CTH  (McGlaun  et  al.  1990)  is  used  at  the  U.S.  Army  Research 
Laboratory  (ARL)  to  model  ballistic  impact  and  penetration  experiments.  The  Johnson-Cook 
strength  model  (Johnson  and  Cook  1983)  is  one  of  several  strength  models  available  in  CTH.  It 
is  an  empirical  model  that  computes  material  flow  stress  as  a  function  of  strain  (work)  hardening, 
strain-rate  hardening,  and  thermal  softening.  The  Johnson-Cook  model  takes  the  following  form: 


^  +  ~8“l(l  +  Clns*)( 

1-r") 

[at 

(1) 


* 


where  A,  B,  C,  m,  and  n  are  constants,  8  is  the  equivalent  plastic  strain,  8*  is  the  strain  rate 
nondimensionalized  by  the  reference  strain  rate  of  1/s,  and  T*  is  the  nondimensional 
temperature.  Parameter  A,  the  initial  (8  »  0)  yield  strength  of  the  material  at  a  plastic  strain  rate 
of  8  =  1/s  and  room  temperature  (298  K),  is  modified  by  a  strain-hardening  factor  (containing 
parameters  B  and  n),  a  strain-rate-hardening  :^tor  (containing  parameter  C),  and  a 
thermal-softening  factor  (containing  parameter  m). 

T*  is  defined  by 


T  = 


T-T, 
T„  -T,  ’ 


(2) 


where  Tr  is  room  temperature  and  Tn,  is  the  melting  temperature  of  the  material,  1,783  K  for 
RHA.  Equation  (2)  is  the  form  used  in  CTH  and  is  valid  for  T,  £  T  ^  Tm,  the  region  of  interest  in 
most  ballistic  apphcations. 


2 


CTH  originally  contained  a  single  set  of  parameters  that  had  been  typically  used  in 
simulations  for  any  thickness  of  RHA.  These  parameters  were  taken  from  one  of  two  data  fits 
for  RHA  presented  in  Gray  et  al.  (1994).  Both  of  these  fits  (which  will  be  discussed)  were 
determined  using  2-in-thick  RHA  that  conformed  to  MIL-A-12560H.  The  fits  resulted  in 
overprediction  of  the  quasi-static  yield  strength  (A  in  equation  [1]).  Their  approach  to 
optimization  was  to  consider  all  parameters  simultaneously.  This  approach  to  fitting  the  data 
resulted  in  a  model  for  the  RHA  that  underpredicted  the  depth  of  penetration  of  several 
experiments;  this  is  discussed  in  more  detail  later.  In  the  present  work,  Johnson-Cook 
parameters  are  developed  for  that  particular  batch  of  2-in-thick  RHA.  The  approach  taken  here 
is  to  fix  the  value  of  A  based  on  the  quasi-static  test  data.  An  optimum  fit  to  the  data  for  each  of 
the  remaining  parameters  is  then  found  individually,  as  suggested  by  Johnson  and  Cook  (1983). 

2.  Dynamic  Data 

Gray  et  al.  (1994)  generated  compressive  stress-strain  data  for  a  variety  of  metals  over  a 
range  of  temperatures  and  strain  rates  using  the  split-Hopkinson  pressure  bar.  The  digital  data 
consisted  of  the  results  of  dynamic  tests  of  2-in-thick  RHA.  At  room  temperature,  tests  were 
conducted  at  four  strain  rates  (0.001;  0.1;  3,500;  and  7,000/s).  At  elevated  temperatures  (473  and 
673  K)  tests  were  conducted  at  a  strain  rate  of  3,000/s.  Strains  were  recorded  firom  near  zero  up 
to  about  0.20. 

To  expedite  processing  time  and  utilize  all  of  the  available  data,  the  digital  data  was  not  used 
directly  to  obtain  the  Johnson-Cook  parameters,  rather  it  was  fit  to  analytical  frmctions  that  were 
suitable  to  the  software  available  for  use  dtuing  this  study.  The  fits  of  the  six  data  sets  are  shown 
graphically  in  Figure  2  and  algebraically  in  Table  1 .  The  frmctions  in  Table  1  are  fits  to  the  RHA 
strength  data  finm  Gray  et  al.  (1994)  and  are  used  to  determine  the  Johnson-Cook  parameters  in 
the  following  analyses.  For  clarity,  yield  strength  predicted  by  die  Johnson-Cook  model  is 
denoted  by  Y  (in  GPa),  whereas  y  (in  GPa)  represents  the  data  fits. 
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Figure  2.  Dynamic  Strength  of  2-in  RHA. 

Table  1.  Dynamic  Strength  Data  for  2-in  RHA 


(3) 


3.  Quasi-Static  Data 

Benck  (1976)  determined  several  properties  for  three  thicknesses  of  RHA.  He  measured  the 
quasi-static  tensile  yield  strength  in  the  three  principal  plate  directions  (in  the  rolling  direction, 
across  the  rolling  direction,  and  through  the  thickness)  at  a  strain  rate  of  0.0003/s.  For  present 
purposes,  these  values  were  averaged  to  obtain  a  representative  isotropic  value.  The 
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compressive  yield  strength  of  the  material  is  then  assumed  to  be  equal  to  the  tensile  yield 
strength;  this  is  only  approximately  true  for  RHA.  The  data  are  presented  in  Table  2  and  include 
unpubhshed  data  for  3/16-in  (4.76  mm)  RHA  (Bruchey  1997). 

The  values  finm  Table  2  and  an  analytical  fit  to  these  data  are  plotted  in  Figure  3.  A 
logarithmic  form  was  chosen;  the  computed  fit  of  the  data  is 

y  =  (-  0. 1428)lnt  +  0.8772 ,  (4) 

where  y  is  the  yield  strength  in  GPa  and  t  is  the  plate  thickness  in  inches. 


Table  2.  Quasi-Static  Yield  Strength  of  RHA 


Plate  Thickness 
(in  [mm]) 

Yield  Strength 
(GPa) 

0.1875 

[4.76] 

1.14 

0.5 

[12.7] 

0.94 

1.5 

_[38.1L 

0.82 

4.0 

[101.6] 

0.69 

4.  Fitting  the  Parameters 


4.1  Parameter  A.  Parameter  A  is  the  yield  strength  at  room  temperature  and  a  strain  rate  of 
1/s.  Equations  (3a)  through  (3d)  in  Table  1  were  interpolated  to  generate  a  fimction  describing 
the  behavior  of  the  2-in  RHA  at  a  strain  rate  of  1/s.  The  resulting  function  is 


y  =  1.5384e°““. 


(5) 


A  comparison  of  equations  (5)  and  (3a)  (Table  1)  shows  a  difference  of  less  than  2%  between  the 
yield  strength  at  e  =  1/s  and  e  =  0.001/s  at  a  strain  of  0.01.  Furthermore,  Benck  and  Robitaille 
(1977)  report  a  difference  of  about  1.1%  for  38-mm  RHA  plate  and  about  0.6%  for  100-mm 
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Figure  3.  Quasi-Static  Yield  Strength  of  RHA. 


RHA  plate  between  the  quasi-static  yield  strength  at  0.0003/s  and  at  0.42/s.  For  the  present 
work,  the  value  of  A  is  approximated  by  the  quasi-static  data  (Table  2  and  equation  [4]).  A  value 
of  A  =  0.78  GPa  for  2-in-thick  RHA  is  obtained  from  equation  (4). 

4.2  Parameters  B  and  n.  The  remaining  parameters  from  equation  (1)  (B,  C,  m,  and  n)  are 
fit  to  the  functions  of  Table  1  by  a  least-squares  technique.  To  fit  the  parameters  B  and  n,  write 
the  first  two  terms  of  equation  (1)  as 


Y  =  A|^l  +  ^e“^ 


(6) 


Equation  (6)  represents  the  yield  strength  at  room  temperature  and  strain  rate  of  1/s,  conditions 
that  render  the  last  two  terms  in  equation  (1)  equal  to  unity.  Equation  (6)  is  rearranged  to 
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Y-A  =  Be“. 


(7) 


Let 


<p  =  ln(Y-A)  (8) 

so  that 

(p  =  nlne  +  b,  (9) 


where  b  =  In  B . 

The  data  at  these  conditions  are  given  by  equation  (5),  which  is  of  the  form 

y  =  Pe“.  (10) 

The  data  must  be  in  the  same  form  as  equation  (8),  so  A  is  subtracted  from  both  sides  of 
equation  (10),  leading  to  the  following  representation  of  the  data: 

^  =  ln(y-A),  (11) 

and 

'T  =  ln(pe“ -a).  (12) 

The  error  incurred  by  ^proximating  the  data  (equation  [11])  with  the  model  (equation  [8])  at  a 
strain  Sj  is  <Pj  -  'Fj .  Subscript  i  represents  an  arbitrary  discretization  of  the  data  into  seven 
strains  covering  the  range  of  the  data  (s  =  0.01, 0.02,  0.04,  0.08, 0.12, 0.16,  and  0.20).  This  was 
done  to  simplify  the  fitting  procedure.  In  the  least-squares  method,  the  error  is  squared  (to  avoid 
having  positive  and  negative  errors  combining  arithmetically  to  reduce  the  total  error)  and 
summed  over  the  range  of  the  data;  the  sum  of  the  squared  errors  is  to  be  minimized: 
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7 

^  ( (pj  -  'Pj  )^  =  minimum.  (13) 

i=l 

The  sum  can  be  minimized  with  respect  to  paramet^  B  and  n  if  the  derivatives  are  set  equal  to 
zero.  That  is. 


(14) 


and 


-T,)’ =0. 

on 


(15) 


Equation  (9)  is  substituted,  and  the  differentiation  is  carried  out,  resulting  in  the  following  two 
equations; 


STilngiSlnej  -lyjSCbSj)^ 
(Ilnsif -7Z(lnEi)^ 


(16) 


and 


Elne, 


(17) 


where  ='i'i(8i)  is  known  (equation  [12]),  and  the  summation  indexes  have  been  omitted  for 

clarity.  The  results  are  B  =  0.78  GPa  and  n  =  0.106.  These  results  minimize  the  error;  this  was 
verified  by  determining  that  the  derivatives  of  equations  (14)  and  (15)  were  positive. 

4.3  Parameter  C.  The  first  two  factors  in  equation  (1)  are 
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(18) 


Y(Ei,E*=l,  T*  =  o)=Afl  +— E”l  =  Si, 

V  A  ; 


where,  for  simplicity,  this  contribution  to  the  strength  is  termed  Sj.  Thus,  for  room  temperature, 
equation  (1)  becomes 

Yj  =Si(l  +  Chis*)(l-0)=Si  +SiClnE‘.  (19) 

To  obtain  a  corresponding  expression  for  the  data,  equations  (3a)  through  (3d)  (Table  1)  are 
used  to  generate  curves  of  stress  vs.  Ins*  for  various  constant  strains.  To  generate  the  curves,  the 
first  of  the  seven  discrete  strains  was  substituted  into  each  of  equations  (3a)  through  (3d) 
(Table  1)  to  generate  stresses  for  each  of  the  four  strain  rates.  The  resulting  stress  was  plotted 
against  Ine* ,  and  the  process  repeated  for  each  of  the  remaining  six  strains,  resulting  in  seven 
curves.  Analytical  expressions  were  fit  to  the  curves;  the  results  are  detailed  in  Table  3. 


Table  3.  Strain-Rate  Dependence  of  Room-Temperatnre  Data 


i 

Strain 

y  =  Ui  +  Vjlne* 

Equation  No. 

Ui 

Vi 

1 

0.01 

1.2587 

0.0035 

(20a) 

2 

0.02 

1.2972 

0.0039 

(20b) 

3 

0.04 

1.3370 

0.0044 

(20c) 

4 

0.08 

1.3780 

0.0050 

(20d) 

5 

0.12 

1.4026 

0.0053 

(20e) 

6 

0.16 

1.4203 

0.0055 

(20f) 

7 

0.20 

1.4342 

0.0057 

(20g) 

These  data  can  be  represented  by  the  form 


yi=S|+S|k|liK*, 


(21) 


where  the  subscript  i  denotes  each  of  the  seven  discrete  strains  in  Table  3.  Equations  (19)  and 
(21)  are  substituted  into  the  least-squares  function 

^i:(Y,-y,f  =0,  (22) 

resulting  in  the  following  equation  for  C: 


C  = 


Isfk, 
Ssf  ■ 


(23) 


This  was  verified  to  be  a  minimum  by  examining  the  derivative  with  respect  to  C  of 
equation  (22).  The  solution  of  equation  (23),  using  all  seven  sets  of  constants  in  Table  3,  is 
C  =  0.0035.  Using  only  the  constants  for  e  =  0.20,  an  upper  limit  (limited  by  the  extent  of  the 
data)  of  C  for  large  strains  is  C  =  0.0040.  The  lower  limit  (using  only  e  =  0.01  constants)  is 
C  =  0.0028. 


4.4  Parameter  m.  Parameter  m  is  determined  by  a  technique  based  on  a  method  described 
by  Johnson  and  Cook  (1983).  Equation  (1)  is  rewritten  in  the  form 

Y,=Y,(l-r”),  (24) 

where  Yj  represents  the  strength  at  temperature  T,  and  the  room-temperature  strength  is 

(25) 


Y,  =A 


l-K^8"j(l-HClnE*). 


The  thermal-softening  factor  is  obtained  from  equation  (24): 
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(26) 


(p  =  lii 


1-- 


■r  J 


SO  that,  also 


<p  =  mlnT*. 


(27) 


Note  that  9  =  cp  (e,  e,  T* ),  but  strain  rate  will  be  held  constant  for  this  analysis.  The  data  must  be 
reduced  to  a  similar  form: 


'T  =  hi 


^  V  "1 

l_Il. 

V  yJ 


(28) 


Only  data  at  a  strain  rate  of  3,000/s  were  used  for  this  analysis.  Since  the  original  data  (Table  1) 
do  not  include  room-temperature  data  at  this  strain  rate,  the  available  room  temperature  data  (at 
strain  rates  of  0.001/s,  0.1/s,  3,500/s,  and  7,000/s)  were  interpolated  to  obtain  values  of  strength 
(yr)  at  a  strain  rate  of  3,000/s. 

Data  in  the  form  yi/yr  were  plotted  vs.  T* ,  and  analytical  fits  were  determined.  Here,  a  good 
fit  to  the  data  was  quadratic: 


^  =  aiT*"+biT*+Ci.  (29) 

Yr 


The  results  of  the  fitting  are  given  in  Table  4. 

Reference  data  is  available  at  293  K  (i.e.,  yJ.  Elevated  temperature  data  is  available  at 

473  K  and  673  K.  Hence,  the  least  squares  optimization  is  fit  for  m  at  two  known  values  of  T* , 
denoted  by  the  subscript). 
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Table  4.  Quadratic  Fit  of  tbe  Temperature  Data 


i 

Strain 

f(T‘)=yT/y,=aiT*'+b,.T*+Ci 

Equation  No. 

ai 

bi 

Ci 

1 

0.01 

“0.2436 

-0.7662 

1.0093 

(30a) 

2 

0.02 

-0.2030 

-0.8029 

1.0057 

(30b) 

3 

0.04 

-0.1627 

-0.8394 

1.0021 

(30c) 

4 

0.08 

-0.1229 

-0.8755 

0.9985 

r  (30d) 

5 

0.12 

-0.0998 

-0.8965 

0.9965 

(30ej 

6 

0.16 

-0.0835 

-0.9114 

0.9950 

(30f) 

7 

0.20 

-0.0708 

-0.9228 

0.9939 

(30g) 

To  obtain  a  globally  (i.e.,  over  the  range  of  temperatures)  optimized  value  for  m,  the 
optimization  formulation  is 


2  7 

SSfej  (31) 

j-i  i=i 

Equation  27  is  substituted,  and  the  result  is  differentiated  with  respect  to  m  and  set  equal  to 
zero.  This  yields  the  following  equation  for  m.  The  summation  indices  have  been  omitted  for 
clarity,  and  'F-  =  'Ey  (sj ,  Tj* )  is  obtained  from  equations  (28)  to  (30). 


(32) 


The  result  of  this  optimization  over  the  range  of  strains  is  m  =  1.07.  Plots  of  the  quadratic  fits  in 
Table  4  show  that  the  curves  approach  a  straight  line  as  s  approaches  a  value  slightly  higher  than 
0.20  (i.e.,  slightly  beyond  the  range  of  the  data).  Hence,  for  strains  slightly  greater  than  0.20, 
m  =  1.  For  small  strains,  equation  (31)  was  solved  for  only  the  8  =  0.01  term  (i  =  1),  with  the 
result  m  =  1.18.  A  value  of  m  =  1.00  was  chosen  for  the  present  work  (see  section  5)  in  order  to 
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favor  the  higher  strains.  This  resulted  in  a  variation  of  less  than  1%  in  depth  of  penetration  when 
compared  to  a  similar  simulation  using  a  value  of  m  =  1.07. 

5.  Numerical  Simulations 

The  continuum  mechanics  code  CTH  was  used  to  run  a  series  of  three-dimensional  (3-D) 
computations  simulating  experimental  ballistic  impacts.  These  simulations  were  run  in  order  to 
validate  the  aforementioned  method  for  determining  Johnson-Cook  fit  parameters  for  RHA.  All 
computations  were  conducted  on  a  Cray  J932  computer  at  ARL’s  Major  Shared  Resource  Center 
(MSRC). 

S.l  Setup.  The  simulations  were  modeled  after  experiments  that  were  perfomied  by 
Enderlein  (1991)  in  which  20-mm-diameter  by  100-mm-long  (1/d  =  5)  X21C  tungsten  alloy  rods 
(right  circular  cylinders)  were  fired  into  stacks  of  four  2.5-in-thick  (63.5-mm)  RHA  plates  at  0° 
obliquity.  These  experiments  were  chosen  as  model  candidates  because  the  2.5-in-thick  RHA 
plates  that  were  used  closely  resembled  the  2-in  (51 -mm)  RHA  for  which  dynamic  strength  data 
were  available,  to  the  extent  that  they  are  in  the  same  thickness  group  in  MIL-A-12560H 
(Figure  1).  Penetrator  impact  velocities  used  in  the  simulation  matrix  corresponded  to  measured 
experimental  striking  velocities  of  1,399  m/s  and  1,616  m/s. 

Figure  4  shows  the  setup  that  was  used  for  each  simulation.  Each  problem  was  modeled  in 
three  dimensions.  The  penetrator  was  inserted  as  a  right  circular  cylinder  with  dimensions  that 
matched  the  experiment.  Eleven  Lagrangian  tracers  were  placed  along  its  centerline  to  monitor 
several  variables,  including  depth  of  penetration.  The  penetrator  was  fired  in  the  x  direction  at 
the  stack  of  RHA  plates.  The  2.5-in  RHA  plates  were  inserted  with  a  0.5-mm  gap  between  each 
other  to  ensure  (numerically)  the  existence  of  an  interface.  The  lateral  dimensions  of  the  RHA 
plates  were  chosen  to  be  8  in  x  8  in  (203  mm  x  203  mm)  in  the  simulations.  Symmetry  was  used 
about  the  XZ  plane  in  order  to  reduce  problem  size.  A  constant  cell  size  of  2  mm  was  selected  to 
finely  resolve  each  problem,  and  this  resulted  in  a  volume  of  approximately  one  million  cells  per 
problem. 
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Figure  4.  Simulation  Setup. 


The  Mie-Griinicsen  equation  of  state  and  the  Johnson-Cook  strength  model  were  used  for  all 
materials.  Default  Johnson-Cook  parameters  were  taken  for  the  penetrator  in  each  simulation. 
These  parameters  are  given  by  Johnson  and  Cook  (1983)  for  90%  tungsten  alloy,  Johnson-Cook 
parameters  were  then  varied  for  the  RHA  in  each  run. 

Six  simulations  were  run  for  each  of  the  two  different  striking  velocities.  Each  simulation  set 
corresponds  to  a  particular  group  of  Johnson-Cook  parameters,  as  shown  in  Table  5.  Set  I 
contains  the  default  parameters  given  in  CTH.  They  are  one  of  the  sets  derived  in  Gray  et  al. 
(1994).  A  nondimensional  temperature  definition  different  from  that  of  equation  (2)  was  used  by 
Gray  et  al,  (1994),  so  using  these  parameters  in  CTH  may  not  be  optimum.  Set  2  is  the  set  of 
parameters  derived  above  for  2-in  RHA.  Sets  3  through  6  are  vaiiations  included  for  information 
and  lend  support  to  the  validity  of  Set  2  for  2 -in  RHA.  These  sets  are  discussed  in  the  following 
sections. 

5.2  Results.  Table  6  compares  the  measured  penetrations  for  the  experiment  and  the  six  sets 
of  simulations.  The  experimental  data  consist  of  a  single  shot  at  each  velocity.  The  depth  of 
penetration  for  each  simulation  was  determined  by  using  the  tracer  at  the  nose  of  the  penetrator. 
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Table  5.  RHA  Parameter  Sets  Evaluated 


Parameter 

Set  V 

Set  2*’ 

Set  3= 

Set  4'* 

Set  5' 

Set  6' 

A  (GPa) 

1.832 

0.78 

0.74 

1.225 

0.9 

0.78 

B  (GPa) 

1.685 

0.78 

0.78 

1.575 

1.305 

0.78 

n 

0.754 

0.106 

0.106 

0.768 

0.90 

0.106 

c 

0.00435 

0.004 

0.004 

0.0049 

0.0575 

0.0891 

m 

0.80 

1.00 

1.00 

1.09 

1.075 

1.00 

“  CTH  default. 

2-in  RHA  (recommended  parameters  for  2-in  RHA). 

®  2.5-m  RHA  (estimated). 

''  Fit  from  Gray  et  al.  (1994)  using  all  strain  rate  data. 

'  Fit  from  Gray  et  al.  (1994)  using  only  high  strain  rate  data. 

^  Current  fitting  technique  using  only  high  strain  rate  data  to  fit  C, 


Table  6.  RHA  Penetrations  in  Millimeters  for  the  Parameter  Sets  Evaluated 


Velocity 

Experiment 

Setl 

Set  2 

Set  3 

Set  4  1 

Set  5 

Set  6 

1,399  m/s 

101 

72 

97 

99 

85 

85 

71 

1,616  m/s 

133 

91 

117 

119 

105 

104 

88 

This  is  illustrated  in  Figure  5,  which  shows  a  time  history  of  the  x  position  of  the  nose  and  tail 
tracers  for  the  l,616--ni/s  2-in  RHA  run. 

Set  1  in  Table  6  shows  that  the  default  Johnson-Cook  RHA  parameters  in  CTFI  lead  to 
underprediction  of  the  depth  of  penetration  for  these  experiments.  In  fact,  the  penetration  is 
underpredicted  by  more  than  28%  for  both  velocities.  There  is  a  significant  improvement  in  the 
penetration  results  when  the  parameters  for  the  fit  to  the  2-in  dynamic  data  are  used  (Set  2  of 
Table  6).  The  penetrations  increase  to  within  4%  of  experiment  for  the  1,399-m/s  case  and  to 
within  12%  of  experiment  for  the  1,616-m/s  case. 

The  underprediction  of  penetration  depth  when  using  the  parameters  for  the  2-in  RHA  was 
expected  because  the  experiment  used  the  softer  2.5-in  RHA  (see  Figure  3).  Had  the  parameters 
for  2.5-in  RHA  been  available,  the  simulation  with  these  parameters  may  have  shown 
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Figure  5.  Nose  and  Tall  Tracer  Histories  for  the  1,616  m/s,  Set  2  Case. 

penetrations  even  closer  to  expCE-iment.  To  demonstrate  this,  the  A  parameter  was  changed  to 
0.74  GPa  to  approximate  the  quasi-static  yield  strength  of  2.5-in  RHA.  As  expected,  the  results 
From  this  simulation  (Set  3)  show  a  slight  improvement  in  predicted  penetration  for  both 
velocities.  For  future  analysis,  dynamic  data  should  be  obtained  for  2.5-in  RHA,  fit  using  the 
aforementioned  method,  and  modeled  in  CTH. 

Set  4  of  Table  6  shows  the  result  of  using  the  Gray  et  al.  (1994)  fit  for  2-in- thick  RHA  that 
considers  the  same  temperature  function  as  CTH,  that  of  equation  (2),  This  set  of  parameters 
predicts  penetration  into  semi-infinite  RHA  better  than  when  using  the  CTH  default  set  (Set  1), 
but  not  as  well  as  when  using  the  parameter  set  derived  here  (Set  2).  Sets  5  and  6  are  discussed 
in  the  next  section. 
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5.3  Discussion.  The  extent  of  the  experimental  variation  may  be  inferred  from  three  shots 
fired  by  Enderlein  (1991)  into  stacks  of  2. 5 -in  RHA  at  30°  obliquity.  The  velocities  of  these 
shots  were  nominally  1,606  m/s  (±2  ra/s),  and  penetrations  varied  by  9%.  In  addition,  a  30° 
oblique  shot  at  1,652  m/s  penetrated  132  mm,  indicating  that  the  data  in  Table  6  for  the 
high-velocity  case  may  be  above  average.  This  experimental  variation  may  reduce  the 
discrepancy  between  the  predicted  and  experimental  values  in  Table  6  for  the  high-velocity  case. 

Several  factors  other  than  RHA  strength  will  affect  depth  of  penetration.  One  possible 
explanation  for  the  slightly  low  CTH  predictions  is  that  the  model  for  the  penetrator  does  not 
precisely  mimic  the  X21C  (93%)  tungsten  alloy  used  in  the  experiments.  The  validity  of  the 
model  was  not  tested  here;  it  is  a  widely  used  model  and  was  chosen  to  facilitate  comparisons  to 
other  work. 

Another  possibility  is  that  the  simulated  interfaces  between  the  plates  do  not  accurately 
imitate  the  experimental  ones  (for  one  reason,  they  are  made  up  of  mixed  cells).  This  may  be  a 
more  important  factor  for  the  high-velocity  cases  because  the  penetrators  get  further  into  the 
second  plate  and  the  model  must  contend  with  a  second  interface  (Figure  6). 

Hardness  of  RHA  (and  therefore  strength)  is  known  to  vary  through  the  thickness  of  the 
plate,  especially  for  the  thicker  stock.  Since  the  RHA  was  modeled  as  homogeneous  for  a  given 
thickness,  this  inhomogeneity  could  lead  to  differences  between  the  experiments  and  the 
simulations. 

Failure  (of  both  the  RHA  and  tungsten)  was  modeled  with  a  basic  threshold  fracture  model 
only  (the  “pfrac”  parameter  in  CTH),  This  parameter  is  typically  chosen  as  the  spall  strength  of 
the  material  (2.5  GPa  for  RHA  and  3.0  GPa  for  tungsten).  This  model  was  kept  constant  during 
this  study  in  order  to  simplify  interpretation  of  results.  No  other  damage  or  fracture  models  were 
used.  However,  fracture  and  damage  are  important  in  ballistic  problems,  as  shown,  for  example, 
by  Raftenberg  (1997a).  This  area  warrants  further  study. 
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Figure  6.  Set  2:  1,616-ni/s  Penetration  Plot. 


CTH  plots  showed  that  the  strain  rate  at  the  penetration  interface  is  as  high  as  lO‘'/s  at  initial 
impact  and  sustains  values  on  the  order  of  [0‘*/s  during  steady-state  penetration,  Raltenberg 
(1997b)  pointed  out  that  the  Johnson-Cook  form  of  strain-rate  dependence  may  not  fit  the 
response  of  RHA  well  lor  this  problem;  C  may  deviate  from  constant  at  strain- rates  greater  than 
lOVs.  This  is  indeed  the  case  with  the  dynamic  data  studied  here  (Table  1).  In  order  to 
investigate  the  effect  of  this  deviation,  two  additional  simulation  sets  were  conducted  (Sets  5  and 
6  in  Tables  5  and  6).  Set  5  uses  the  Gray  et  al.  (1994)  fit  of  only  the  high  strain-rate  data.  The 
result  was  virtually  unchanged  from  Set  4,  with  the  Gray  et  al.  (1994)  fit  using  all  of  the  data.  A 
similar  attempt  was  made  here,  wherein  the  parameter  C  was  fit  using  only  the  high  strain-rate 
data  (3,500/s  and  7,000/s);  a  value  of  0.0891  resulted.  All  other  parameters  were  unchanged 
from  Set  2.  The  resulting  penetration  depths  at  both  velocities  were  severely  underpredicted. 
This  underprediction  is  not  surprising,  since  making  such  a  change  to  C,  which  is  a  constant,  will 
strengthen  the  RHA  at  all  strain  rates  throughout  the  block.  The  least-squares  method  used  for 
the  strength  data  in  Table  5  (Set  2)  selected  an  optimum  value  for  C  over  the  full  range  of  strain 
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rates.  This  should  provide  a  reasonable  value  for  a  variety  of  applications,  but  further  study  of  C 
is  warranted. 

6.  Conclusion 

A  method  for  fitting  static  and  dynamic  strength  data  for  RHA  in  order  to  obtain  parameters 
for  the  Johnson-Cook  constitutive  model  has  been  presented.  This  method  can  be  used  to  obtain 
parameters  for  available  thicknesses  of  RHA,  provided  dynamic  data  exist.  Parameters  were 
obtained  from  2-m  RHA  data,  tested  using  CTH,  and  shown  to  significantly  improve  predictions 
for  a  particular  set  of  experimental  data  when  compared  to  predictions  obtained  using  two  sets  of 
previously  published  parameters,  including  the  CTH  default  set. 
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Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 
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